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ABSTRACT 

The electrical power subsystem is the heart of a satellite; it provide sa permanent power source to all subsystems 
in the satellite. The power sources are either primary source or secondary source, primary power source is the 
main source that supplies power to satellite subsystems in sunlight period, and secondary power source is a 
back-up source that supplies the load during eclipse periods.This paper focuses on using multi -junction solar 
array as primary power source and the lithium ion capacitor as secondary power source. Size of solar array 
and lithium ion capacitor for micro-satellites power systemare presented. A case study and its simulation using 
Matlab/Simulink are analyzed. The results show that lithium-ion capacitor can be used as an energy storage 
alone or without need of chemical battery beside. 
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Abbreviation: 



SA Solar Array 

Lie lithium-ion capacitor 

LIB lithium-ionbattery 

LEO Low Earth Orbit 

EDLC Electric Double Layer Capacitor 

P^onsumed daily average power consumed 

P^enerated daily avcragc power generated 

BCS Body Coordinate System 

OCS Local Orbital Coordinate System 

RAAN Right Ascension of Ascending Node 

LAN Longitude of Ascending Node angle 

ESR Equivalent series resistance 



Nomenclature: 



P(t) 


Generated power at any instant is given 


So 


Irradiance coefficient 


riT 


Overall efficiency of power subsystem 


A 


SA area 


n" 


SA normal vector 


S" 


Sun vector direction 


a 


Sun angle; between " and ^ 


risA 


SA efficiency of conversion 


iIpss 


Power subsystem efficiency 


Kff 


S A fill factor 


Ktemp 


Temperature coefficient 


Koeg 


Degradation coefficient 


e 


Elevation angle; between the Z-axis and " 


o 


Azimuth angle; between the X-axis and " 


V 


nu angle; between ( ^ ), and the normal to the orbital plane (Y-axis) 


P 


Satellite in-orbit position angle; between ^ and the noon line. 




Inclination angle; between the orbital plane and the equatorial plane 


i(t) 


i at any instant 


io 


Initial i at moment of launching 




Rate of change i 




LAN angle; between X-axis, and the line of intersection between the orbital plane and the equatorial plane 


Q(t) 


The LAN angle at any instant 


Oo 


Initial LAN at moment of launching 
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Rate of change of LAN 

Longitude angle of the sun (RAAN); between X-axis, and the Hne of intersection between the ecliptic plane and the equatorial 
plane 

The RAAN angle at any instant 
Initial RAAN at moment of launching 



^(t) 




P 

S 




I 

dt 



dv 



H(t) 
Ho 
H' 
5 



G 



S 



Rate of change of RAAN 

Inclination of the ecliptic plane to the equatorial plane (23.4523o ) 
The attitude of the orbit at any instant 
Initial altitude at moment of launching 
Rate of change of altitude 

Roll angle; rotation of the satellite around the X-axis 

Pitch angle; the rotation of the satellite around the Y-axis 

Yaw angle; the rotation of the satellite around the Z-axis 

Change in voltage from the Vmax to Vmin 

Maximum working voltage 

Minimum working voltage 

Current withdraw by capacitor 

The discharge time corresponding to drawn current 

Total capacitance of the system 

Total resistance 

Voltage drop caused by the ESR 
Cell capacity 
ESR of cell 

Total numbers of parallel cells 
Total numbers of series cells. 



I. 



INTRODUCTION 



A satellite consists of two main parts. The first one is called payload and the second is platform [1]. 
Payload is the satellite main purpose, andplatform contains communication, telemetry, on-board computer, 
electrical power and storage system and attitude and orbit control. Micro -Satellites have shown great capability 
in recent years in the field of remote sensing, such as earth imaging, disaster monitoring, military application, 
and deep space exploration. The small satellites have the abilityof re-visiting times for remote sensing missions 
and providing larger instantaneous coverage areas. The main problem in using the micro -satellite is the power 
generated not big enough to supply any device with high power for short duration like thrusters. In Low Earth 
Orbit (LEO) spacecraft where secondary source are typically cycled with a 100-min cycle periods consisting of a 
65-min charge and a 35-min discharge [2]. 

Solar Array (SA) is the most widely and common source of electrical power in space. SA convert sunlight 
into electricity for Earth orbiting satellites. Power requirements in tens of watts to several kilowatts, and a life 
time from a few months to 15 years can be met with a SA. The output power of SA depends on temperature, 
space environment, irradiance, dirt, shadow, and so on. It produces energy according to its operational 
conditions, location, and environmental conditions [1]. 

Recently, the multi-junction solar array is used because it has higher efficiency and performance than 
silicon cell. It generates higher power compared to silicon cell, and has smaller temperature coefficient so it can 
also be used in interplanetary missions. [3,4]. 

Chemical batteries is the most realistic choice for energy storage device, because its performance is 
well-known, trusted, and used in many and many space missions. Now, the Lithium ion Batteries (LIBs) is the 
common secondary power source in satellites as its cell voltage is about 3.6 vwhich is much higher than NiCd 
and NiH2 by three times, have a good temperature range,have a long cycle life, do not suffer from memory 
effect and have high specific energy above 130 Wh/Kg.[5]The disadvantage of using the LIB are it is strong 
dependent on the temperature and its performance deteriorated at low temperature so it have a negative impact 
on the thermal system, also is in order to keep the batteries functioning for the entire mission duration which 
contain more than 30,000 charging-discharge cycles, the Depth of Discharge (DoD) must be worked at low 
value about 30% as number of cycles of battery is inversely proportional with the DoD, so the power subsystem 
will use more number of cells to compensate the decrease in DoD. So this problem should be solved. 

Today, super-capacitors(SCs) have become a good alternative to chemical batteries as 
rechargeable energy storage medium. Unlike batteries, super-capacitors provide very huge numbers of cycles 
(can reach to 1,000,000 cycles). In addition, the high specific power nature of super-capacitors will enable 
increase the payload and other subsystems power on board such as actuators and thrusters. Some additional 
benefits are also wide operating temperature range, and can easily measure the state of charge. 

Super-Capacitors or ultra-capacitors also called electrochemical double layer capacitors are energy 
storage devices that store electric energy in the double layer between a high surface area electrode and an 
electrolyte [6]. As compared to batteries, there is no chemical reactions between the electrodes (not all SCs types) 
that is a purely physical phenomenon. Therefore, electrical energy storage produced in absence of any chemical 
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reactions [7]. Due to the purely physical operation, the SC scan provide a large amount of cycling without any 
problem. Furthermore, the performance of SC with decreasing temperature does not deteriorate strongly as 
compared to the LIB [6]. SCs offer several advantages over secondary batteries, such as high power capability, 
and high charge-discharge efficiency, etc. [8]. However their low specific energy (range of 5 Wh/kg for cell) [9] 
limits their applications as secondary energy storage source. Here SCs are used to complement main energy 
sources, such as Solar array (SA), and batteries, when there is an instantaneous demand for relatively large 
power. 

A lithium-ion capacitor (LIC) is hybrid capacitor that combines the features of SCs and lithium-ion batteries 
(LIBs), thus offering relatively large specific energy without losing the advantages of conventional SCs[10], and 
give higher discharge period. 

Super-Capacitors are an energy storage material that have a very huge capacitance in hundreds and thousands 
farads, and high specific power [11]. 

They can be classified into three major types; Electric Double Layer Capacitor (EDLC), pseudo -capacitance and 
hybrid super-capacitors [12], as follows: 

A. Electric Double Layer Capacitor (EDLC) 

The energy is stored in the double layer formed at the electrode-electrolyte interface. Activated carbon is 
mainly used as the electrode material because it has low cost, and high surface area. EDLC is used in vehicles 
applications and it is the most common typethat has high specific power and low specific energy. 

B. Pseudo-capacitor 

Thiscapacitor uses oxidation-reduction reaction and stores the chargefaradaically by transferring of electron 
chargebetween the electrode and the electrolyte. A pseudo-capacitor has a chemical reaction at the electrode. 
Ituses metal oxide or conducting polymer as electrode materials. Its advantagesare having greater capacitance, 
and much higher specific energy than the EDLC. This technology is mainly in the research phase and is not 
commercially available [12]. 

C. Hybrid Capacitor 

This capacitor has characteristicsof combination the features of LIB and EDLCs. The positive electrode has 
activated carbon and the negative electrode has carbon material that is doped with lithium ion. LIC is the most 
popular hybrid capacitors. The advantage is having much larger specific energy, the cell voltage is higher, and 
have large discharged time period. Sothe overall performance of LIC is better than EDLC [10]. 



II. SOLAR ARRAY SIZING 

The sizing of SA means calculating the SA area needed that covers the power consumption. 
The SA generated power at any instant is given by the following equation: 



P(t) = Sq X iij X A X cosa(t) 



Where ^- is irradiance coefficient, is the overall efficiency, - is the SA area, - is the sun angle 

between the normal to SA vector (^) and the direction of sun vector ^S)^ "si the SA efficiency of 

conversion from light energy to electrical energy , 'jP-S power subsystem efficiency, F- is the SA fill factor, 

^te _p is the Temperature coefficient, and ^0 - g is the degradation coefficient. In the design, the Power Supply 
Subsystem (PSS) must meet the requirement; where the daily average power generated is greater than the daily 

average power consumed at any moment during the lifetime PGane rated ^Consunisd [13, 14]. 
Moreover, if the SA area is required for a given consumed power, we will re-arrange equation (1). 

P 

A = ^ — 

So X iiT X COS! a,i (3) 

If the ^- P, and^- remain unchanged, the SA area is affected only by sun angle. Therefore, the problem here is 
to find the value of sun angle and compute the SA sizing. 

Since the dot product gives the angle between two vectors, therefore we will use it to find the sun angle value. So 
^ and S is given by the following equation: 
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n S = \n\\S\ cos(a) (4) 

Assume that the norm 1^1 and IS| are a unit vectors, and re-arrange the previous equation. 

cos(g:) = n ■ S (5) 

First, we have to calculate the values ^ and ^ then compute the sun angle. 
Calculation of ^ : 

^ will be calculated in Body Coordinate System (BCS) where selection of system axes is arbitrary according to 
the mission. 

According to Fig. 1, ^ described in BCS: 

Secs = sin(e)cos(<|>.)i -J- 3in(e) s[n(<|>)j + cos(9)k 





(a) 



(b) 



In matrix form 



Figure 1 (a) SA at an arbitrary position and (b) in BCS 



"bcs — 



siniej cosiC©:) 
sin(e)sin(0) 
cos(e) 



(7) 



Where - is the elevation angle between the Z-axis and ^ , and ^ is the azimuth angle between the X-axis 
and n. 

Calculation of ^ : 

^ will calculate in Local Orbital Coordinate System (OCS) where selection of system axes are; the origin is 
located at the center of mass of the satellite, X-axis is the positive direction of orbital flight, Y-axis is the normal 
to orbital plane, and Z-axis is the direction from the satellite to the Earth[15] 

According to Fig. (2), S described in OCS: 

Sqcs = ~ siii(v)s]n(p) I + cos(v) S]ii(v)cos(p)k (g) 

In matrix form 

— siniv J sin(p) [I 
Socs = cosCv.) J 

sin(v)co3(p) J Ik J (^9^ 

Where - is angle between ( S ), and the normal to the orbital plane (Y-axis), and - is the satellite in-orbit 
position anglebetween ^ and the noon line. 

The importance of v is that, this angle shows us the worst-case day over the satellite lifetime. The v angle is not 
to be conftised with the sun angle (a), often used to define the sunlight incidence angle on the solar array, which 
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could be canted to catch the sun normally in case v is not 90°. The power generation is proportional to cos (a), 
so when ais 0° the power generation is maximum and 90° gives zero power. 
The V angle is given by the following equation: 



v(t} = -Ts{u^\ 



sini(t ) sinQ(t) cosi|/(t) 
+ sine cosiCt) 3iin|r(t) 
^+ c 0 s E si n i (t) c 0 3 flCt.) sin t) / ( i q) 



Where iis inclination angle of the orbital plane to the equatorial plane, ^ is longitude angle of ascending node 
(LAN) 




Figure 2 (a) OCS, (b) ^ in OCS 

Fig. 3, ^" is Longitude angle of the sun in ecliptic plane "Right Ascension of Ascending Node" (RAAN), and ^ 
is inclination angle of the ecliptic plane to the equatorial plane Fig. 4. 

The inclination angle (i) at any instant, The LAN angle ( ^) at any instant. The RAAN angle ( - ) at any 
instant, the attitude of the orbit (H) at any instantare given by the following equations respectively : 



i(t) = io + xt 
H(t) = Ho + H' xt 



(11) 
(12) 
(13) 
(14) 



Where ^ a is inclination angle at moment of launching, and ^ ' is rate of change inclination, is the initial LAN 
at moment of launching, is rate of change of LAN, '^d is the initial RAAN at moment of launching, and ^ * 
rate change of RAAN, -^a is the initial altitude at moment of launching, and ^ is the rate of change of altitude. 



orbital plane 




Equatorial plane 

Ascending Node 
Figure 3 Orbital plane of Satellite 
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z 




OCS. The rotation matrices are given by the following equations: 



Rv = 



10 0 

0 cos(6) — sin(S) 
0 sini(6) cos(6]) 

cos(6) 0 sinT(6]) 
0 1 ' 0 

— sini'ej 0 cos(6) 

cos(ji) sm(n) 0 
sini'i-i) C03(|i) 0 
0 0 1 



(15) 



(16) 



(17) 

■Iocs ~ I^Z^^Y^^X^BCS 



(18) 



Where ^ is theroll angle; rotation of the satellite around the X-axis, is pitch angle; the rotation of the 
satellite around the Y-axis, and .^^ is yaw angle; the rotation of the satellite around the Z-axis. The 
multiplication order of matrices is important. 

CO SI. a ' is called an illumination coefficient (Kni), not used in the calculation directly but it need some 
modification e.g. if > 90°, then cos (^^ ) will be a negative value and since either power or area cannot be a 
negative values, therefore the working region of sun angle from 0° to 90° otherwise cos( ) will be zero. 



III. 



SUPER-CAPACITOR SIZING 



The SC voltage contain two components; a capacitive component represents the voltage change due to 
the change in energy within the SC, and a resistive component represents the voltage change due to the 
equivalent series resistance (ESR). 

The capacitive component, the resistive component, and the SC total voltage components aregiven by the 
following equations respectively: 



d V = I X 



dt 



(19) 



V = IxR, 



(20) 
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(21) 



Wheredv is change in voltage from the maximum voltage (Vmax) to minimum voltage (Vmin), I is the current, dt 
is the discharge time corresponding to withdraw current, CTotai is the total capacitance of the system, V is voltage 
drop caused by the internal resistance and RTotai is the total resistance. 

The capacity of the module (CTotai). and the resistance of the module (RxotaOaregiven by the following equations 
respectively: 

P 

(22) 

5 

(23) 

Where Cceiiis cell capacity, Rceiiis the equivalent series resistance per cell, P is the number of parallel cells, and 
S is the number of series cells. 

To determine numbers cells are required in series, divide the maximum applied voltage by the cell voltage. 
To determine numbers cells are required in parallel, we will re -arrange the equation (20) to be: 



P = I X 

dv(^ + RcEii) 



(24) 



So, now we know the all parameters in sizing a bank of SC.There are different methods in sizing SC depending 
on using the previous equations either by using the energy in SC or by using the SC time constant or by ignoring 
the effect ofESR. 

A. Energy in SC 

This method works well for low power application where the loss due to ESR is minimal. Since the energy 
storage in SC is equal to ¥2 CxotaiV^, therefore the total capacitance is: 

CJotal = 2f/C:T^mai^2 - 7j^min^2 ) ) (25) 

Where E is energy of SC needed can be calculated from multiple the average power by working time. 

B. SC time constant 

We will assume SC time constant is equal 1 second or if we use the specific SC then, we will know the Rceii and 
Cceii- Since RxC=l, R=l/C, then use it in equation (20) to find Cxotai- 

C. Ignore the ESR effect 

ESR to find the value of Cxotai then use the equation (20) as a check for the design after choosing the proper SC 
cell. The total capacitance calculated by this method is 

dt 

^rotc; — I ~ 

-1-" (26) 

Note that the difference between these methods are small, the second, and third method did not take into account 
any losses due to ESR, we typically recommend oversizing the cell by adding a cell in series and parallel but not 
necessary. 
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IV. 



CASE STUDY 



Now we will present the result in our case study .The case study based on using the Matlab, Mathcad, 
and STK software programs to simulate the behavior of power subsystem in satellite in the desired orbit. The 
SA is fixed around Y-axis,total numbers of panels are eight as shown in Fig. 5. 




Jody Z 



Figure 5 Satellite in the desired orbit 



The results are based on some acceptable assumptions; the earth has spherical shape, the sun light is parallel and 
this lead to lack of penumbras, during study one orbit the sun, the earth, and the satellite assumed to be fixed 
during the whole revolution period, and the orbit has circular shape. 

V. RESULTS 

As shown in_Fig. (3) the satellite has daily average power consumed (Pconsumed) equal 81.58 w , and max power 
equal 233.8 w in duration 119.4 sec 
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Figure 6 Load profile during a typical day 
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300 




Time (min) 
Figure 8 Generated power in one orbit 
Solid line is average power generated, dotted line is average power consumed 
The orbital inputs parameters shown in table (1) in appendix 

A. SA sizing 

The SA parameters is shown in table (2), and the SA data sheet in appendix (1). 

Using equation (7) to find ^EIS , equations (11-14) to find v angle, equation (9) to find'^O^S , then equation 

(18) to find *^0IS , at this step we know the value of values ^ and S so apply equation (5) to find the value of 
cos (a) and its correction factor as mentioned before it will 0.31. In equation(3), then area needed to generated 
82W is equal to 1.36 m^ then use the safety factor 10% then the daily average power generated (Pgenerated) equal 
90.20 w which is larger than Pconsumedby 10% and the max power equal 283.38 w which is larger than the max 
power consumed by 17.4% as shown in Fig. (7) 

B. SC sizing 

As shown in Fig. (7) or (8) the shadow duration is equal to half of total orbit duration. The SC parameters is 
shown in table (3), and the SA data sheet in appendix (2). 

We will use the Vceii at 3.5 V, so S = 9, in equation (23) P =9 and add another one in case of cell failure, so P 
=10 T= n +T V2 T =PA^ T =PA^ T=2 94 A dv=V -V =4 V 

1 ) Energy in SC 

In equation (25), CTotai=2376.21 F, calculate Cceii— 

2138.59F, so we will choose Cceii— 2300 F so calculate Ctom 
again it will be 2555. 56F then make a check using equation (21) the dv new =3.4V <dv old ,so the design is 
accepted, the energy storage in SC is equal to 79.5 IWh. 

2 ) SC time constant 

According to table (3), RC=1.38, in equation (21), CTotaF2172.28 F, CceiFl955.05 F, then choose 
Cceii=2300 F it will give the same result as above. 

3) Ignore the ESR effect 

In equation (25), CTotaF2171.27 F, CceiFl954.14 F, then choose CceiF2300 F it will give the same result as 
above. 

Note that at the beginning in design, there are difference between methods of design but at the end, the results 
are same. Now will check the whether this design can withstand during max current in max power duration. 
Imax=Pmax/Vbus=8.36 A, dt =119.4 scc, therefore the dv=0.4 V, so the design is good. 

VI. Appendix 



Table 1 Orbital inputs j 


parameters 


Parameter 


Value 


Launching day 


21/3/2014 


Ho 


680 Km 


Maximum change in altitude 


OKm 


Satellite lifetime (Tiife) 


5 years 


io 


98.1115° 


Maximum change in inclination 


0° 


Tlan 


10.5 AM 


elevation angle,0 


0° 


azimuth angle, O 


0° 
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Table 2SA parameters 



Parameters 


Value 


-^consumed 


82 W 


So 


1360 W/m^ 


SA efficiency, rjsA 


29.5% 


SA fill factor, KpF 


85% 


SA efficiency of PSS, rjpss 


90% 


Temperature coefficient, KTemp 


80% 


Degradation coefficient, Koeg 


80% 



Table 3SC parameters 



Parameters 


Value 


-^consumed 


82 W 


T shadow 


49.175 min 


V 

max 


30 V 


V 

min 


26 V 


^bus 


28 V 


Vceii(rated) 


3.8 V 


Ccell 


2300 F 


ESRcell 


0.6 mQ 



Table 4 Solar array Emcore,ZTJ Photovoltaic cell, space photovoltaics 



Parameter 


Value 


BOL Efficiency at Maximum 


29.5% 


Power Point 




Voc 


680 Km 


Jsc 


OKm 




5 years 


Jmp 


98.1115° 



Table 5Lithium ion capacitor Ultimo lithium ion capacitor (Prismatic series) 



Parameter 


Value 


Vceii(max.) 


3.8 V 


Vceii(min.) 


2.2 V 


Ccell 


2300 F 


ESRcell 


0.6 mQ 


Specific energy 


8 Wh/kg 



VII. CONCLUSION 

In this work, a power supply consists of solar array and Lithium-ion capacitor is designed to electrify 
micro-satellites. Solar Array is used to supply the loads and charging the storage source. Lithium-ion capacitor 
supplies the load in case of eclipse periods. A complete design for the SA is suggested based on the total 
consumed load power. A check is done to ensure that the generated power is more than the consumed power 
under any case whether in average or in maximum and during the eclipse duration. Super capacitor is designed 
based on the previous information. The results show that the LIC can supply the load properly by storing the 
sufficient energy during the day. The SCs in space application is look promising as it has huge number of cycles, 
wide temperature range, high efficiency, and its problem of specific energy can be overcome. 
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